function A=fill_function(Phi)
                     
A=cell(4,4);

A(:,1)=Phi(:,1);
A(:,2)=Phi(:,2);
A(:,3)=Phi(:,3);
A(:,4)=Phi(:,4);
end